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Abstract. We solve the time-dependent dynamics of ax- 
isymmetric, general relativistic MHD winds from rotating 
neutron stars. The mass loss rate is obtained self consis- 
tently as a solution of the MHD equations, subject to a fi- 
nite thermal pressure at the stellar surface. Conditions are 
chosen to be representative of the neutrino driven phase 
in newly born magnetars, which have been considered as 
a possible engine for GRBs. We compute the angular mo- 
mentum and energy losses as a function of a and compare 
them with the analytic expectation from the classical the- 
ory of pulsar winds. We observe the convergence to the 
force-free limit in the energy loss and we study the evolu- 
tion of the closed zone for increasing magnetization. Re- 
sults also show that the dipolar magnetic field and the 
presence of a closed zone do not modify significantly the 
acceleration and collimation properties of the wind. 



1. Introduction 

Magnetically dominated outflows from stars and accre- 
tion discs, are a key element in understanding the evo- 
lution and properties of these objects. Magnetic outflows 
are ubiquitous in powering a variety of astrophysical sys- 
tems, by converting the rotational energy into wind ki- 
netic energy. Rotationally powered pulsars are among the 
best example of relativistic outflows produced by a fast 
compact rotator. The high magnetic field at the surface 
(B ~ 10 10 " 13 G) and short rotation period (P - 0.001 - 1 
s) cause particles injection. If the charge density in the 
magnctosphere of the neutron star, due to pair creation, 
exceeds the Goldreich- Julian value then conditions can be 
modeled in terms of ideal RMHD. 

Observations of the interaction of the pulsar wind with 
the surrounding SNR, tell us that far away from the source 
the wind is ultrarelativisitic (7 ~ 10 4-6 ) and weakly mag- 
netized (jArons 2004[) , in contrast with conditions inside 
the Light Cylinder, where 7 ~ 100, and magnetization is 
high. This implies that the high Lorentz factor is achieved 
by almost completely converting magnetic energy into ki- 
netic energy 



In this sense, pulsars constitute the prototype of rela- 
tivistic accelerators: a plasma is assumed at the base with 
small outflow velocity (compared to the asymptotic veloc- 
ity), embedded in a strong magnetic field. Magnetorota- 
tion then accelerates the plasma and converts the mag- 
netic energy into relativistic motion. 

In the case of pulsars, the plasma density at the sur- 
face of the neutron star, drops to very low values, and the 
magnetosphere is strongly magnetically dominated. Such 
conditions are close to force-free, and the equations that 
describe the flow dynamics become stiff, thus making im- 
possible for present day codes to model such environment, 
in terms of MHD. Due to the presence of gaps at the foot 
points of open field lines, plasma is injected in the mag- 
netosphere at a supersonic speed, quite differently from 
what happen in normal stars, where injection is subsonic. 

However low magnetized neutron star winds are of 
interest primarily in understanding the physics of very 
young neutron stars. In the seconds after collapse and ex- 
plosion, the neutron star is hot (surface thermal speed 
~ 0.1 — 0.03c ), and has an "extended atmosphere". 
This proto-ncutron star cools and contracts on its Kelvin- 
Helmholtz timescale (tkh ~ 10 — 100 s), radiating its grav- 
itational binding energy (~10 53 ergs) in neutrinos of all 
species (jPons et al. 1999|) . The cooling epoch is accompa- 
nied by a thermal wind, driven by neutrino energy deposi- 
tion, which emerges into the post-supernova-shock ejecta 
( [Thompson et al. 200l] >. 

After birth, the thermal pressure at the edge of the pro- 
toneutron star surface decreases sharply as the neutrino 
luminosity decreases on a timescale tkh , as the star cools 
and deleptonizes. At some point during the cooling epoch 
the magnetic energy density will exceed the thermal pres- 
sure, and the wind region will become increasingly mag- 
netically dominated. For magnetar-strcngth surface fields 
(Bo ~ 10 14 — 10 15 G) the magnetic field dominates the 
wind dynamics from just a few seconds after the super- 
nova dThompson 20031 ). 

If magnetars are born with millisecond rotation peri- 



ods (Thompson & Duncan 1993), their rotational energy 
can be very large (~2 x 10 52 ergs) relative to the en- 
ergy of the supernova explosion (~10 51 ergs) and be- 
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Fig. 1. Wind structure for a monopolar magnetic field, a = 200, fl = 0.14. Left: poloidal magnetic field lines, and 
ratio \B,p/B r \; Alfvenic (dashed) and fast (dotted) surfaces. Center: Lorentz factor, ~ sin(0). Right: energy flux in 
sin 2 (9). From Bucciantini et al.(2006) 
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cause the spindown timescale can be short for larger value 
of the Alfvenic radius, these proto-magnetar winds have 
been considered as a mechanism for producing hyper- 
energetic supernovae ( Thompson et al. 2004[ ). Because the 
wind becomes increasingly magnetically-dominated, the 
outflow must eventually become relativistic. For this rea- 
son, proto-magnetar winds have also been considered as 
a possible central engine for long-duration GRBs (e.g., 



Thompson et al. 2004 ) 



In the simple ID radial model for relativistic winds 



from compact rotators by Michel (1969) solutions can 
be parameterized in terms of a = Q 2 & 2 /M, where ft 
is the rotation rate, <f> the magnetic flux and M the 
mass flux, a is the maximum Lorentz factor the wind 
can achieve if all magnetic energy is converted into ki- 
netic energy, and at the base of a pulsar it is ~ 10 6 . 



The results by Michel (1969) show that the asymptotic 
Lorentz factor tends to a L/i . This is known as the tr — 7 
problem (|Arons 200 4). If the flow is allowed to diverge 
more than radially then acceleration is more efficient 
(Begclm an fc Li 1994| ). The presence of a magnetic field 
can cause collimation of the flow along the rotation axis, 
and divergence along the equator, so in principle one might 
expect magnetic collimation to be also responsible for the 
acceleration of the wind. However when the flow speed 
reaches values close to c the Coulomb term in the elec- 
tromagnetic force cannot be neglected and it balances the 
magnetic hoop stress. Solutions of the relativistic Grad- 
Shafranov equation for a monopolar field in the trans-fast 
regime (Beskin et al. 1998) have shown that acceleration 
is too slow and the asymptotic results of the Michel model 
apply also in 2D axisymmetric geometry, meaning that the 
flow remains radial. 

However, no analytic multidimensional solution of the 
full relativistic MHD equations exists that extends from 
the surface of the neutron star to large distances. Am- 
bigiuties remain about injection conditions, which in the 
case of a newly-born hot magnetar, depend non-trivially 
on the temperature, gravity, and rotation. Even simplified 



models consider only supersonic injection, but none focus 
on the possibility of a thermal wind at the base. Moreover 
the use of a monopolar field might prove not accurate, 
given than a more realistic dipole should be assumed in 
the close magnetophere. 



2. Wind acceleration 

The first axisymmetric numerical simulations of the 
acceleration of a relativisitc wind were presented by 
where 



Bogovalov (1999) where the case of a zero pressure 
plasma in a monopolar magnetic field, with uniform in- 
jection conditions at the base of the star was investigated. 
Results show that as soon as 7 > 5 collimation and accel- 
eration of the wind become marginal. 

Our recent investigation (Bucciantini et al. 2006) has 
focused on the study of the acceleration of an axisym- 
metric relativisitc wind in the case the pressure is small 
compared with the rest mass energy density, but not neg- 
ligible. If acceleration is efficient, the wind can reaches 
high Lorentz factors (7 ~ 100 — 300) and kinetic luminos- 
ity (~ 10 50 erg s _1 ), perhaps enough to trigger a GRB. 
The relativistic MHD equations, together with Maxwell 
equations, are solved using the shock-capturing code for 



relativistic MHD developed by Del Zanna et al. (2003) 



The code has been modified to solve the equations in the 
Schwarzschild metric, in order to properly account the ef- 
fects of gravity. In order to model the heating and cooling 
processes using an ideal gas equation of state we adopt an 
adiabatic coefficient L = 1.1 (almost isothermal wind). A 
hot atmosphere is assumed at the base of the star with 
p/ pc 2 ~ 0.01. Density, pressure, rotation rate and the ra- 
dial component of the magnetic field are fixed at the base 
of the neutron star, which is assumed to have a canoni- 
cal radius Rns = 10 km. The system is let to relax to 
a steady state configuration. We verify that it takes less 
than 5 rotations for the mass flux and energy losses to 
converge. 
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Fig. 2. Wind structure for a dipolar magnetic field, a = 20, O = 0.14. Left: poloidal magnetic field lines, and ratio 
\B,p/B r \; Alfvenic (dashed) and slow (dotted) surfaces. Center: same as left, note the monopolar structure at larger 
radii. Fast (dotted) surface. Right: Lorentz factor, ~ sin {&). From Bucciantini et al.(2006) 



In analogy with the ID model, solutions in the 2D 
case can be parametrized by using the average value of 
a over all stream lines. Results in the case of a monopo- 
lar magnetic field show that for a > 10 the flow is close 
to radial, the terminal 7 is < a, pointing to inefficient 
acceleration. Moreover the energy and angular momen- 
tum fluxes rapidly approach the force- free values. Even 
if v << c inside the Light Cylinder, asymptotic results 



agree with the flow structure found by Bogovalov (1999) 
Energy flux scales as sin 2 (0), and 7 oc sin(0), in agree 
mcnt with the solution of the exact monopole. Mass flux 
is also shaped by rotation. On the equator it appears to 
be 5 times higher than in the non rotating case. On the 
contrary due to collimation in the polar region the mass 
loss rate along the axis is about twice smaller than in the 
non rotating case. Such variation is due to the change in 
centrifugal support at different latitudes, and decreases as 
the period increses. Results suggest that a rapidly rotat- 
ing magnctar will not be able to simulataneusly produce a 
collimated and relativisitc outflow. Relativisitic velocities 
can be reached only for a » 1, and in this case most of 
the energy would be concentrated in the equatorial plane 
instead that in a jet. Simulations have been extended to 
large distances (~ 10000-Rats) showing that for a > 10 
most of the energy is in the equatorial plane. We observe 
the formation of a collimated narrow jet, that is unre- 
solved in our grid, but such structure appears not to be 
energetically significant. 

We computed solutions also in the more realistic cases 
of dipolar magnetic field at the surface of the NS. Due 
to stability issues in numerical RMHD it was not possi- 
ble to study cases with a > 20. In the monopolar case, 
no current sheet is present. In the case of a dipole cur- 
rent sheets are present both at the edge of the closed zone 
and at the equator, and numerical resistivity is important. 
In particular in the equatorial plane, where the mangetic 
field vanishes, reconnection (due to numerical resistivity) 
takes place and plasmoids are formed. These plasmoids 
extend about 5-10 degrees around the equator, and cause 
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Fig. 3. Energy and angular momentum losses for a 
monopolar (upper panel) and dipolar (lower panel) mag- 
netic field as a function of a a (a averaged on open field 
lines). Squares indicate cases with different rotation rate. 
Dashed-lines represent the low-cr analytic theory, dotted 
lines are the force-free limits. The solid line is a power law 
fit for the monopolar field. When normalized on open field 
lines, the results of the dipole follow the same curve of the 



monopole. From Bucciantini et al.(2006) 
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fluctuation of about 10-20% in the losses. Even if such 
values depend on numerical resistivity, this suggests that 
reconnection at the Y-point might manifest itself observa- 
tionally. 

To obtain steady solutions we suppressed resistivity at 
the equator, thus enhancing the instability of the code. 
Moreover, in the case of a dipole, the magnetic and flow 
surfaces are not aligned with the grid, and in this case in- 
version from conservative to primitive varibles is not stable 
for B 2 /pc 2 > 100. 

Despite the presence of a closed zone, and the fact 
that within the Light Cylinder the field structure is close 
to a dipole, we found that solutions can again be par- 
materized in term of cr, now defined as an average done 
only on open streamlines. In this case we find that energy 
and angular momentum losses scale as a monopolc with 
the same a. Again as soon as cr > 1 solutions for the en- 
ergy and angular momentum losses approach the force- free 
limit dContopoulos et al. 1999 ), outside the Light Cylin- 
der, and for a > 10 the asymptotic structure resembles 
the exact monopole, with energy flux peaking at the equa- 
tor and a higher equatorial Lorentz factor. Similar results 



are found also by Komissarov (2006) We are able to re- 
lated the amount of open magnetic flux to the distance 
of the Y-point from the surface of the star. Extrapolat- 
ing to the force-free limits, our result is in agreement with 



Gruzinov (2005) 



Contrary to expectations dContopoulos et al. 1999D we 
also find that the location of the Y-point is inside the Light 
Cylinder. This is due to the presence of a finite pressure 
plasma at the base of the star. In this case the extent of 
the closed zone is limited by the pressure equilibrium at 
the Y-point. We also verified that the ratio of the Y-point 
radius to the Light Cylinder Ry /Rl drops with decreasing 
il, as one would expenct for a braking index less than 
3. This implies that even if for large a the asymptotic 
wind structure is close to the force- free monopole solution, 
the spindown time can be shorter, and the total energy 
losses larger, that what is inferred by applying the dipole 
formula. 



3. Conclusions 

Recent numerical results have confirmed theoretical pre- 
dictions about the acceleration of relativisitc winds from 
neutron stars. The conversion of magnetic to kinetic en- 
ergy appears to be minimal, the flow closely resembles the 
expectation from the exact monopole solution, there is no 
evidence for the formation of a collimated energetic jet. 
The use of a more realistic dipolar field does not seem 
to change significanlty the wind properties far from the 
Light Cylinder. On the contrary the energy and angular 
momentum losses seem to depend only on the integrated 
value of cr, and, for a > 1, on the amount of open magnetic 
flux. 



Regarding the possibility that a magnetar could trig- 
ger a GRB, we conclude that, in ideal MHD, acceleration 
to relativistic speeds requires very high cr, much higher 
than what is thought to characterized newly-born magne- 
tars. Moreover there is no evidence of energy collimation 
at high magnetization. Even if dissipation of the current 
sheet could be invoked to accelerate the wind, it is not 
obvious that a jet will form, it is more likely that the en- 
ergy will still be concentrated in the equatorial plane, as 
is suggested by X-ray images of Pulsar Wind Nebulae. 

Even if pulsar winds are not thermally driven at the 
base, they are cold, and the injection conditions depend 
on the physics of gaps, our results have several speculative 
implications for rotation powered pulsars, where a is high. 
The fact that the Y-point is interior to the light cylinder 
(and the Alfven radius), and that Ry/Rl decreases as 
decreases, suggests that observed braking indices less than 
3 might be a consequence of pressure equilibrium. Even if 
globally, thermal pressure is not likely to be relevant for 
classical pulsars, magnetic dissipation at and near the Y 
point might cause pressure forces to be significant in this 
localized region. 

The appearance of outwardly propagating plasmoids 
at and beyond the Y-point as a consequence of (numerical) 
magnetic dissipation raises the intriguing possibility that 
noise in pulsar spindown might arise from instability of 
the magnetospheric currents. Torque fluctuations might 
give rise to noise and limit cycles in the observed phase of 
pulsars' subpulses. 
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